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We show how the analytical approach of Zakharov and 
Rubenchik [Sov. Phys. JETP 38, 494 (1974)] to modu- 
lational instability (MI) of solitary waves in the nonlinear 
Schroedinger equation (NLS) can be generalised for models 
with two phase symmetries. MI of three-wave parametric 
spatial solitons due to group velocity dispersion (GVD) is 
investigated as a typical example of such models. We reveal 
a new branch of neck instability, which dominates the usual 
snake type MI found for normal GVD. The resultant nonlin- 
ear evolution is thereby qualitatively different from cases with 
only a single phase symmetry. 



One of the central issues of solitary wave theory is 
the question of stability ||l|,|| . Stability of solitary waves 
('solitons') in nonintegrable Hamiltonian models is often 
governed by the derivative of a certain invariant with re- 
spect to an associated parameter of the solution [Cl-pl. 
For instance, positivity of the derivative of the total 
energy with respect to the nonlinearity-induced wave- 
number shift is sufficient for stability of bright solitary 
waves in such fundamental processes as self-action of ra- 
diation in media with intensity dependent refractive in- 
dex and in degenerate three-wave mixing in quadratic 
nonlinear media H. 

A wide range of parametric processes conserve not only 
the total energy of the interacting waves but also other 
energy invariants. These can be, for instance, energies of 
individual waves or energy unbalance of two waves. Typ- 
ical examples of such processes are non-degenerate three- 
wave P-u| and four- wave mixings ||]. By Noether's the- 
orem every integral of the motion in Hamiltonian models 
originates from a corresponding symmetry property. En- 
ergy invariants are often associated with phase (gauge) 
symmetries. Although many issues of the dynamics of the 
multi-wave solitons in models with several symmetries 
still remain to be understood, their stability threshold in 
many situations is given by a zero of the determinant of 
the Jacobi matrix H^K^QiH [H,||J7|- Here Qi and Hj are 
respectively the integrals of motion and their associated 
parameters. 

Soliton dynamics in a space where the soliton is lo- 
calised in some dimensions but extended in one or more 
others, e.g. consideration of a soliton stripe, raises the 
problem of modulational instability (MI) along the extra 
dimensions ||l| . MI of solitary waves has been studied in 



mariy fields, including plasma physics |I|, fluid dynam- 
ics |l],pX]| and optics [p],pT|-[l^ . For instance in optics the 
above-mentioned extra dimensions might be associated 
with diffraction and/or group velocity dispersion (GVD). 
An analytical approach to the low-frequency limit of MI 
for bright solitary waves was originally developed by Za- 
kharov and Rubenchik Q for the generalised NLS equa- 
tion, which has a single phase symmetry. This approach 
is based on asymptotic expansion near neutrally stable 
eigenmodes (Goldstone modes) of the solitary wave ex- 
citation spectrum. The presence of such neutral modes 
is directly linked to the symmetries of the model. It fol- 
lows that systems with a single phase symmetry should 
be qualitatively similar to the NLS case p], where the 
soliton always shows MI in the extra dimension. The 
unstable mode is of neck type (the soliton stripe breaks 
into a chain of spots) or snake type (the stripe distorts in 
zig-zag fashion) depending on the relative signs of the dis- 
persive terms in the localised and extended dimensions. 
All previous studies of solitary wave MI |]l],|0-15| have 



been restricted to situations with a single phase symme- 
try, and all do indeed exhibit NLS-like MI. For example, 
in degenerate three- wave mixing (3WM), which has a 
single phase symmetry, the first analytical results on MI 
of solitary waves reported by Kanashov and Rubenchik 
|ll| and recently extended and supported by numerical 
results 12,11], typically show the neck/snake scenario. 

The influence of extra phase symmetries on MI of soli- 
tary waves is still an open issue and is the main sub- 
ject of the present Letter. We will concentrate on non- 
degenerate 3WM as a typical and practically important 
example of a solitonic model with two phase symmetries, 
motivated by recent theoretical ]^,0] and experimental 
P,n advances in the study of quadratic optical solitons. 

In this Letter we show that MI of solitary waves in non- 
degenerate 3WM reveals a new branch of neck-type insta- 
bility which can be ascribed to the extra phase symmetry 
and its corresponding neutral mode. We show that in me- 
dia with normal GVD this new instability dominates the 
snake-mode responsible for MI in the corresponding de- 
generate 3WM model. Our theoretical approach and the 
phenomena predicted by it should extend to other soli- 
tonic systems of broad interest with similar symmetry 
properties, e.g. to incoherently coupled NLS equations 
[|l6| and non-degenerate 4WM |||. 

In non-degenerate 3WM the evolution of suitably nor- 
malised slowly varying field envelopes Em {m = 1,2,3) 
of three waves with carrier frequencies LOm (wi + W2 — lo^) 



propagating in a dispersive and diffractive quadratic non- {vi, V2, vs^ 



linear medium can be modeled 
tern of dimensionless equations: 



{wi,W2,W3)'^ and 



by the following sys- 



id.Ei + aidlEi + 71 9? £^1 + £^2 £^3 = 0, 

id,E2 + a2dlE2 + 725(2-^2 + ElE^ = 0, (1) 

id,E:i + a^idlE^i + 73^2^3 + E1E2 ^ fSE^. 

Transverse x, longitudinal z and retarded time t co- 
ordinates are respectively measured in units of a suit- 
able beam width, diffraction length and GVD parame- 
ter. The am and 7™ are diffraction and dispersion co- 
efficients referred to these scales, and the wave- vector 
mismatch is characterised by /3. We neglect spatial and 
temporal walk-off effects, implicitly assuming that either 
their spatial and temporal scales are much longer than 
those associated with diffraction and GVD, or that walk- 
off is compensated by special techniques, as in recent 
experiments on temporal solitons in degenerate 3WM 
[ p7t . Henceforth we make the experimentally appropriate 
choice Q;i,2 = 2a3 = 0.5. We also assume all jm either 
negative (normal GVD) or positive (anomalous GVD), 
leaving the case of mixed dispersion for future work. For 
the sake of simplicity we have restricted our model to one 
transverse dimension, as in a planar waveguide. 

Suppressing the time derivatives for the moment, Eqs. 
(1) have a family of non-diffracting solitonic solutions 
Ejn = A,„(x)e'('''"^+'^'"\ where the Am are real, ^1^2 = 
K ± 6 are positive parameters, n^ — 2k > —j3 and 
4'i,2 = V i V"! '/'3 = 2(/?, where ip and '0 are arbitrary 
real constants. In general functions Am{x) must be found 
numerically or approximated variationally ||6|,[7| . The free 
choice of (/3 and "0 implies two phase (gauge) symmetries, 
which by Noether's theorem leads to two conserved quan- 
tities. These are the total energy Q = Qi -I-Q2 + 2Q3, and 
energy unbalance Qu = Qi — Q2, or equivalent combina- 
tions of the Qm — J dx\Em\'^- Note that the degenerate 
case forces Ei — E2, and thus V' = 0, which suppresses 
one phase symmetry. 

Our primary aim here is to study temporal MI due 
to GVD of these spatial solitons. This is most interest- 
ing and important when they are spatially stable. Their 
stability against purely spatial perturbations has been 
studied P,0|, yielding a stability boundary which in our 
notation is given by dsQd,^Qu — d^^QdsQu- Spatially 
stable domain is in fact almost the entire domain of soli- 
ton existence, excluding only a small range of k, S values 
with /? < 1^,0 . Close to this region there is also a small 
domain of bistability H. Therefore the existence and 
nature of temporal MI is the important question for al- 
most all parameter values, and in particular for the entire 
half-space with /3 > 0. 

To study MI due to GVD we seek solutions of 
Eqs. (1) in the form of spatial solitons weakly mod- 
ulated in time at frequency 17 > 0: Em = (Am(x) + 
(Ki(a;,z) -Kiiy™(x,z))cosm)e*('"'"^+"^'"l Setting Vm = 
Vm^^^ , Wm = Wme^^ , wc obtain two eigenvalue prob- 
lems L-L^v — X^v and L^L^w = X^w, where v = 



L± = 


±Li A3 ±A2 

A3 ±L2 ±Ai 

_ ±^2 ±Ai ±L3 


O^mdl - 7mf^^ 


- ^m, where ^1^2 



Lm = otmdl - 7mri^ - ^„, whcrc ^1^2 = K ± (5 and ^3 = 
2k + (3. Note that L± are self-adjoint, so L_L+ and 
L-^-L- are adjoint operators with identical spectra. It 
is thus enough to consider the spectrum of one of these 
operators, e.g. L+L_. We are particularly interested in 
the discrete spectrum, with eigenfunctions exponentially 
decaying at x —^ ±cxd. 

In general the stability problem can only be solved 
numerically, but for small absolute values of A we can 
obtain some analytical results. Our two phase symme- 
tries, plus the Galilean one, generate three neutrally sta- 
ble (A = 0) eigenmodes of L^L^ at O = 0. These 
are: w^ = (^1,^2, 2^3)^, w^ = (^1,-^2,0)"^, w^ = 
{xAi,xA2,2xA3)'^ . Assuming that 17 <IC I we can ex- 
press w as a linear combination of Wip, w^, Wx and then 
use an asymptotic approach to find eigenvectors and cor- 
responding eigenvalues. This approach can be applied 
only if all other eigenmodes have eigenvalues obeying the 
condition [A| > fi. In practice this only excludes a small 
neighborhood of the spatial stability boundary discussed 
above. 

Three eigenvalue pairs ±A are obtained from solvabil- 
ity conditions of the first order problems. One, associated 
with the asymmetric eigenvector Wx, obeys 
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— — dx y JmidxAm) 



(2) 



Clearly the asymmetric mode is unstable for normal 
GVD, which corresponds to the snake instabilities found 
in NLS §,|l§ and degenerate 3WM |l|,|l| models. 

The other two eigenvalue pairs are associated with 
linear combinations of the spatially symmetric vectors 



a 



ip i^Lp 



C,pw^, and thus with neck- type instabilities. 



They are the roots of 



aX^ + hVt^X^ + cVt'^ = 0, 



(3) 



where a = {dsQdf^Qu ~ dKQdsQu)/2, b = 9kQ(7i(9i + 
72Q2) + dsQuiliQi + 72Q2 + 473Q3) + {d^Qu + 
d5Q){'^2Q2-liQi), andc= -8(7i72(5iQ2+7273Q2Q3 + 
7i73QiQ3)- These expressions are quite complicated, but 
yield some important general results. Clearly c is nega- 
tive when all 7^, of the same sign. Since a > throughout 
the spatially monostable domain, it follows that the two 
roots A^ are always real and of opposite sign, so that 
there is always an unstable neck-type mode. Thus we 
establish coexistence and competition of neck and snake 
instabilities for normal GVD. This is a novel feature of 
the present model, quite different from previous analyt- 
ical results for NLS §,0 and degenerate 3WM |ll|,|^. 



where the snake instabihty is the only one for normal 
GVD. This is because only one symmetric neutral mode 
exists in models with a single phase symmetry, and it 
generates instability only for anomalous GVD. In non- 
degenerate 3WM there are two neck modes, one stable 
for normal GVD and unstable for anomalous GVD, and 
vice versa, since b is odd in the jm- 

Simple analytic expressions for growth rates of these 
neck modes can be obtained in several special cases, e.g. 
in the case of second harmonic generation (wi = uj2)- Set- 
ting 7i — 72 and 6 = (readily achieved in experiment 
g) the two eigenmodes have either C^ = or C^ = 0, 
with eigenvalues 



x2 



27ir2' 



2 Ql 



K 



40^ 



hiQi + 273Q3). (4) 



For 6 = 0, df^Q is positive and dsQi is negative in the 
spatially monostable domain. Thus the novel neck insta- 
bility for 71 < can be directly attributed to the gauge 
symmetry in the differential phase ip and its associated 
neutral mode w^. On the other hand X^p is associated 
with the usual neck MI for anomalous GVD in models 
with a single gauge symmetry &4l31 . The expression for 
A^ holds also (when Q« = 0) for other solitonic models 
with a differential phase symmetry. Note, however, that 
dsQi can generally have either sign, leading to instability 
with either normal or anomalous GVD. 

Solving the eigenvalue problem numerically, we find 
that in low-frequency limit the instability growth rates 
precisely match those predicted by our perturbation the- 
ory, see Fig. 1(a), (b). As fi is increased each MI gain 
curve reaches a maximum and then decreases. A typical 
example of the maximal MI growth rate vs (3 is presented 
in Fig. 1(c). Similar plots for Qu 7^ and across wide 
range of 7,„ values show the same behaviour ||l^] . Thus 
we conclude that for normal GVD the new neck insta- 
bility strongly dominates the snake one. Note that its 
growth rate is maximised, as Fig. 1(a), (b) illustrate, for 
Qu = 0. For normal GVD the unstable eigenfunctions 
become weakly confined and develop oscillating tails as 
fl increases. Because this increases computer demand, we 
have plotted growth rates in Fig. 1(a) only for il values 
corresponding to well-localised eigenmodes. Physically, 
broader eigenmodes have weaker overlap with the soli- 
ton, and hence lesser gain. 

Spatial profiles of the symmetric eigenfunctions at 
maximum gain (f2 = Qmax) are presented in Fig. 2. 
Despite this being well beyond the perturbative limit in 
which expressions (4) apply, the novel neck MI eigen- 
mode still has qualitatively the same form as w-^, i.e. 
wi = —W2, u'3 = 0, indicating that the ijj phase symme- 
try underlies the instability through the whole range of 
51. Similarly, the ustable neck mode for anomalous GVD 
is evidently associated with the (p symmetry. 

To test our linear stability analysis and study the non- 
linear evolution we performed an extensive series of com- 
puter simulations of the system (1) with initial conditions 



in form of a soliton stripe perturbed by spatio-temporal 
white noise of order 1%. Typical simulation results are 
presented in Fig. 3 and they fully support our predic- 
tions. We chose the size of the computational window 
in the time domain to be l87T/fl„iax, and the initial soli- 
ton stripe rapidly develops nine humps, in accord with 
the stability analysis. During further evolution the mod- 
ulated stripe forms into a train of pulses which either 
spread (normal GVD) or form a persistent chain of three- 
wave optical bullets (anomalous GVD) . Due to the initial 
noise, modes from a band of frequencies close to ^max are 
able to grow and compete, and hence the modulations in 
Fig. 3 are somewhat irregular. 

A striking difference between Figs. 3(a,c) is that the 
initially imposed translational symmetry of the solitary 
stripe along the time dimension is broken in different 
ways. For normal GVD interleaved intensity peaks of 
El and E2 are formed, while for anomalous GVD the 
intensity peaks coincide. (Each amplitude is modulated 
with period ~ 2TT/Qmax-) 

This difference is directly related with the spatial form 
of the most unstable eigenvectors. In the case of nor- 
mal GVD wi and W2 are out of phase and wa = 0, 
see Fig. 2(a), leading to the interleaving. Since E1E2 
drives E3, the intensity profile of the second harmonic 
becomes modulated with period n/Qmax, see Fig. 3(a3). 
Because the overlap of the three fields is diminished by 
this evolution, mutual trapping becomes impossible and 
the whole structure eventually spreads through diffrac- 
tion and dispersion, see Fig. 3(b). For anomalous GVD, 
all three components of most unstable eigenvector are in 
phase, see Fig. 2(b), and thus all three intensities be- 
come modulated with the same temporal period, see Fig. 
3 (ci), (C2), (ca). This provides conditions for mutual 
self-trapping of the filaments, see Fig. 3(d). 

The predicted instabilities can be experimentally ob- 
served in a nonlinear material longer than the MI gain 
length Ig ~ A~^ and with pulse width order of several 
271 /rimax or more. Following Ref. |13], typical values for 



KTP Ig - 1cm and 27r/fi^ 



10 s. Fig. 1(c) shows 



that the new neck instability has the highest gain and is 
thus most easily observable. Artificially birefrigent semi- 
conductor materials, which are highly nonlinear, seem 
quite promising, based on a recent experiments on 3WM 
|19[ . Waveguides containing Bragg structures might be 
very suitable for observation of the MI phenomena pre- 
dicted here because of their large and controllable dis- 
persion |20[ . 

In summary, we have analysed and described dispersive 
MI of spatial solitons due to non-degenerate 3WM. Using 
this model as a typical example we generalised a previous 
analytical approach to MI Q to models with two phase 
symmetries. We found that the extra neutral mode asso- 
ciated with the additional phase symmetry gives rise to a 
new branch of MI. This is symmetric (of neck type), and 
is found to dominate the asymmetric (snake) instability 
which is the only MI for normal GVD in systems possess- 
ing just one phase symmetry. This result enables a new 



understanding of the dynamics of multi-component soli- 
tary waves in terms of their phase symmetry properties. 
The MI phenomena which we have described for 3WM 
are likely to be generic in other solitonic and nonlinear 
wave models with two phase symmetries. 
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FIG. 1. (a-b) Instability growth rate vsQ: Q = 65, /3 = 0. 
Thick (thin) lines are for k = 2, 5 = 0, Q„ = (k = 2.075, 
(5 = 1.525, Qu = —36). Full (dashed) lines correspond to neck 
(snake) MI. Dotted lines are perturbative results, (a) Normal 
dispersion: 71,2 = 273 = —0.5. (b) Anomalous dispersion: 
71,2 = 273 = 0.5. (c) MI growth rate at il = fimaa; vs /3 for 
K = 2, (5 = 0. Full (dashed) lines correspond to neck (snake) 



MI for 71,2 = 273 
71,2 = 273 = 0.5 



-0.5; dot-dashed line to neck MI for 
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FIG. 2. Most unstable eigenmodes {Q — flmax) for 
l3 — S — 0, K = 2. (a) Normal dispersion: 71^2 ~ 273 — —0.5; 
(b) Anomalous dispersion: 71,2 ~ 273 — 0.5. 
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FIG. 3. Development of noise-induced instability of spa- 
tial soliton stripe: k = 2, 5 = /3 = 0. Left (right) panels for 
71,2 = 273 = -0.5 (71,2 = 273 = 0.5). (a„) \Em\ at z = 2.7, 
(b) \Ei\ &t z = 4.5, (c„) \Em\ at z = 5.4, (d) |_Bi| at z = 10.8. 



